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^ ! Abstract 

o 

I We study the exact solution of the anisotropic spin-1/2 Heisenberg chain 

i with a boundary magnetic field in the region where the bulk excitations are 

o> ■ 

gapless. It is shown that near the boundary a bound state is created which is 



underscreened when 7 < vr/S and which at low temperatures behaves like a 
single spin weakly coupled to the bulk. The IR fixed point in this case belongs 



j-^ ■ to the universality class of the underscreened anisotropic Kondo model. We 



argue that the same fixed point appears in the boundary sine-Gordon model 
when the scaling dimension of the boundary term 2/3 < A < 1 . 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.20. Dx 

Problems about the influence of impurities and imperfections on the behaviour 
of one dimensional strongly correlated systems are attracting growing attention. 
So far the main interest has been concentrated on effects of potential scattering on 
the behaviour of electrons in Luttinger liquids^'^. It is widely believed that there 
the conductance of a one-dimensional chain vanishes at T = even if only a single 
impurity is present, provided the electron-electron interaction is repulsive ^. The 
best studied model deals with spinless interacting electrons and in the continuous 
limit is equivalent to the so-called boundary sine-Gordon (BSG) model. There 
are strong reasons to believe that this model is related to the problem of the 
spin-1/2 Heisenberg chain with a boundary magnetic field'^: 



H = -y 

2 ^ 

n=l 



\i-(^n(^n+l - ^n+l^n ) + COS 7«+i 



- h2a% (1) 
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We shall discuss this relation below. The model (1), however, is interesting in 
its own right. At hi^2 = it describes an experimentally accessible situation of a 
one-dimensional magnet where some spins are replaced by non-magnetic ions. 

The model (1) is exactly solvable by the Bethe ansatz'^'^. The Bethe ansatz 
equations are 

M 

[ei{Ua)f^ e2sAUa)e2S2iUa) = H ^^{Ua - Ub)e2{Ua + Ub) (2) 



b=l,bjta 



^ sinh[7(^ - in/2)] 
^ ' sinh[7(ii + in/2)] ^ ' 



exp{2i'ySj) = ci 

and the energy is given by 



1 , fhj ^ 
— in — -|- cos 7 

27 J 



(4) 



^ ^ 1 d 1 

where 

- 27rsin7 

u — u 

7 

The quantities 5'i_2 are defined in such a way that at /i = 25" = 77/7— 1 and 
277/7— 1 /i— >oo. Solutions Ua and —Ua describe the same eigenstate. 
Since we are interested in boundary effects we have to compare the free energy 
of the open chain with that of a chain with periodic boundary conditions. In 
order to make this comparison we shall rewrite the Bethe ansatz equations in 
a form where they are maximally similar to the equations for the XXZ model 
with the periodic boundary conditions. Following Ref. 3, we define a new set of 
rapidities Va such that 



Va = 



Ua, a = l,2,...M 

-U2M-a+l, a = M+l,...2M 



Then Eqs.(2) become 



r r \m ( x / ^sinh[27(w„ - i/2)] 

[eiW] --sAva)e2sAva) ^^^^^^^^^ ^ .^^^j - n e.(.. - V,) (6) 

The last term on the left hand side is introduced to compensate the term with 
vi, — —Va now present on the right hand side. This term can be rewritten as 

sinh[27(t;„ - i/2)] 



sinh[27(^;„ + i/2)] 



ei{va)ei-T,/^{va) (7) 
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The term ei{v) can be included into the bulk part. Finally, we get the following 
Bethe ansatz equations: 

2M 
6=1 

with the energy equal to 

1 - 1 d 1 

In these notations the Bethe ansatz equations look similar to the equations for 
the periodic XXZ spin-1/2 Heisenberg chain with an inserted impurity spin^ or 
the equations for an anisotropic Kondo model with the large coupling constant 
(see, for example, Ref. 7). The analogy becomes complete when one of the 
boundaries is free {h2 = 0, 25*2 = 7r/7 — 1). Then the two phases on the left hand 
side of Eq. (j^) cancel and the resulting Bethe ansatz equation is the same as that 
for the Kondo model with impurity spin 5*1. Except for the special case hi oo 
this spin is never equal to 1/2. Since the XXZ Heisenberg chain with an impurity 
spin is very similar to the Kondo problem^, where the singlet ground state is 
guaranteed to exist only if the impurity spin is equal to 1/2, one can expect here 
non-analytic contributions from the boundary. As we shall see later, this does 
occur when < 7 < 7r/2 when the IR fixed point of the model (1) coincides with 
the fixed point of the underscreened Kondo model. At ir/2 < j < n the boundary 
spins are completely screened at T ^ 0. 

Below we shall consider the thermodynamical properties of the model (1). In 
order to simplify the calculations we shall do it at the special points (a) 7 = tt/z/ 
and (b) 7 = 7r(l — z^^^), where the solutions of Eqs.(H) have especially simple 
classification^. The equations for the distribution functions of rapidities coincide 
for the both cases, but the energy has different signs. Since the results for a 
finite boundary magnetic field do not differ qualitatively from those for the open 
boundary conditions and the latter case is more physical, we shall concentrate on 
the problem with open boundary conditions and make only few remarks about 
the general case when it will be appropriate. 

We emphasise that we are interested in boundary effects, but not in finite size 
effects. The latter ones disappear when the length of the chain becomes infinite. 
Neglecting finite size effects we get the following thermodynamic equations for 
the open boundary conditions: 

en = -^Jr]s{v)6n,i + Ts * In (1 + e^"'^/^) (l + e^"+^/^ 
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+T6n,,-2S * In (l + e^«/^) , n = 1, ...v - 2 (10) 



e,_i = Ts * In (l + e^^-^/^) + ^i/ (11) 
eo = Ts * In (l + e^-^'^/^) - (12) 
Ffe^ifc = -2ArT dt;s(t;) In (l + e^^/^) (13) 

^boundary -^1 ~l~ ^2 

Fi = -T dvs{v) In (l + e^-i/^) (l + e^°/^) 

F2 = -T^~dt;s(t;)ln(l + e^^/^) (14) 

where 77 = 1 for 7 < 7r/2 and 77 = — lfor7>7r/2 and H is the magnetic field in 
the bulk. The sign * stands for convolution 

roc 

f*g{v) = / f{v -v')g{v')dv' 



and s{v) = [4 cosh(7rt>/2)]~-^. 

As we have said, one can expect these equations to be similar to the equa- 
tions for the Kondo problem. There is one important difference, however: the 



non-trivial part of the boundary free energy F2 (0) is equal to 1/2 of the free 
energy of the Kondo impurity. This is, of course, due to the restriction that only 
symmetric distribution of rapidities are allowed. Eqs. (10-14) are also similar to 
the thermodynamic equations for BSG problem derived in Refs. (9, 10) from 
the bootstrap solution of Ref. 11. The anisotropy parameter 7 is related to the 
scaling dimension of the boundary term in BSG problem: 

A = l-7/7r (15) 

Here the differences are more important. The least important one is that the 
BSG equations contain the boundary energy scale Tg. This difference by setting 
Tb equal to the ultraviolet cut-off A ^ J. What is more important, however, is 
the fact that the boundary free energy in BSG model has a different form: 

Fbsg = -T r dvs[v - - HX/Tb)] In (l + e^-^'^) (16) 

At i7 = when e^-i = eo and A = Tb these two expressions are equivalent, but 
at 7^ one can expect differences. These differences are due to the assymetry 
between solitons and antisolitons in BSG. 

Since the obtained thermodynamic equations are very similar to those for 
BSG model and the latter were analysed for 77 = —1 in Refs. (9, 10), we shell 
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concentrate on the case = 7 <7r/2. In this case the ground state consists 
of real f 's, i.e. the only non- vanishing energy at T = is ei. 

Analytical solutions are available for asymptotics of e„(a;) at f — 5> 0, oo (see, 
for example, Refs.(7, 8)). At large temperatures T » J the free energy is 
determined by the asymptotics aX v ^ +oo where 

(17) 



(l + e^"/^) (l + e^"-!/^ 
such that we have 



sinh vH/2T 



^boundary ^ 2 



sinh HI2T 
sinh vH/2T 



(18) 



^ In I 

sinh H/2T 

At small temperatures the leading contribution comes from the region i> — > 
where e„ {v ^ 1) are again almost constant and exp(ei/T') is small. Then the 
corrections can be determined from the expansion in exp(ei/T'): 

gr,{v) ^ In (l + e-^'^)/^) = g^^^ + g^:\v) + ... 

= 21n$(n) (n = 2,3,...i/-2), 
g^^l^ = ln[l + ex.^{vH/2T)^{v - 2)], = ln[l + ex.^{-vH /2T)^{v - 2)] 

^" ^^""^ " ^(2)$(n) ^^^"^ + ^^""^ * ^'^""^ ~ ^^'^ ~ * ^'^""^^ ^^^^ 

. . sinhnz/iJ/2(z/ — 1)T , , sinh[(z/ — n)ti;/2l , 
= siuh.g/2(l-l)r - "-f^' = .uhl(.-lW2] <^°) 

Using these expressions we get the following expansion for the free energy: 

Fboundary ^ - - I) - T dvf{v) In (l + 6^^^^)/^) (21) 

where A = 1 — l/u — 1 — 7/77 and 

. _ tanh(co'/2) 
~ sinh[(i/-l)c^/2] 

In order to estimate the second term in Eq.(21) we use the crude approximation: 

In (1 + e^i(^)/^) ^ce(v-^ HJ/T)^ (22) 

where C is a constant. 

The first result is that the ground state has a finite entropy 5'(0) = — | ln(A~-^ — 
1) which corresponds to the the half of the entropy of the underscreened spin 
(5* — 1/2). A careful analysis shows that this entropy disappears at A < 1/2. 
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The ratio of boundary contributions to the partition functions in the ultraviolet 
and the infrared limits is 

ZyylZiR = (23) 

The same expression remains valid at vr > 7 > 7r/2^. 

From the obtained expression for the free energy we derive the following asymp- 
totics for the specific heat and the magnetic susceptibility at u > 2, H 0: 



X - 



24(t/- 1)T 



l + +... (24) 

~ (T/ 7)2/(^^-1) (25) 



where 5 is a constant. 

Now we shall calculate the boundary contribution to the overall magnetic 
moment M^"""'^'^'^?' and the average value of spin on the boundary (Sf) at T = 0. 
At zero temperature the only non- vanishing density is the density of real f 's p{v). 
From the Bethe ansatz equations (8) we derive the integral equation for p{v): 

rB I 

/ duAii{v - u)p{u) = s * Aii{v) + —ui{v) (26) 

J-B 2l\ 

E = -^-JN dvs*Au{v)p{v) (27) 

where 

, , , ^ sinh(z/ - 2)cj/2 , , sinh(z/ - 2S)cu/2 
smh uuj/Z smh uoo/Z 

and the limit B is determined by the magnetic field in the bulk such that at it is 
infinite at if = 0. Let us consider the case h = first. At H << J the Fredholm 
equation ( ^61) can be treated as a Wiener- Hopf equation. The latter can be solved 
giving the leading asymptotics of the magnetization. The result is 

j^toun,ary ^ J_ T ^^ CXp [- ln( J/if )] 1 siuh 

47r J- oo X — iO )(a;) sinh[Ax/(l — A)] coshx 

- ^tan(7r/A):i;ii±l/^(i7/J)2(^-^-i) 

(H/J) (28) 



2v/^ ' ' ^r(l/A-l/2) 

1 r[l + l/2(l- A)] 



v^cos[7r/2(l - A)] r[l + A/2(l - A)] 
where G^'^^x) are functions analytical in the upper (lower) half-plane of x: 

GW(-x) = G(-)(x) = ^ (29) 

r(i + ^)r(i/2) 



6 



At A > 2/3(z/ > 3) the non-analytic contribution becomes dominant. At A = 2/3 
there is a double pole and we have a marg inal situation with M'"^'"^'^y H\nH. 
Comparing Eqs. ( pSj ) and ( p5D we see that H has the same scaling dimension as 
T which is quite natural since if we have an unscreened spin in our theory. 

It is also interesting to know whether the underscreened spins are situated 
directly on the boundaries. To find this out we calculate the average value of the 
boundary spin explicitly at if = 0, T = 0. From Eqs.p^) and ( p7[ ) we find the 
/ii-dependent part of the ground state energy: 

^_h,_ j U^.M d^ sinh[.-2g(M]^ ^^^^ 
2 71 Jo sinh ucu cosh to 

Differentiating this expression in hi and using the definition of S (1) we get 

(5r) = -|f k=o = ^-4/ du;-^ = (31) 

ohi 2 TT'' Jo smhu; 

and 

d'^E 4 r°° duu'^ta.nhLj 

Xboundary or 2 '^1=0 q '■ 7 / '■ w 1 n~ K'^'^) 

oni 7'^sm7JJo smh(7rti;/7) 

Thus we see that the susceptibility of the boundary spins is finite. This means 
that the boundary spins do not participate in the bound state responsible for the 
singularities in the free energy. 

In conclusion, as we have suggested above, at 7 < 7r/3 the boundary contribu- 
tion to the free energy is non-analytic in T and H and the boundary free energy 
is equal to 1/2 of the free energy of the corresponding underscreened Kondo 
model. Since the XXZ chain is similar to BSG model, we expect the same type 
of behaviour for BSG in the area 1/2 < A < 1. 

The authors gratefully acknowledge important conversations with P. Coleman, 
M. Evans, L. loffe, E. Fradkin and I. Kogan. 
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